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Abstract 



We derive the approximate pressure profiles, density profiles, and temperature profiles of an 
atmosphere, also called barometric formulas. Our variant of a derivation goes beyond the 
common standard exercise of a thermodynamics lecture, where commonly the discussion of 
the underlying physical assumptions is missed. We depart from the Navier-Stokes equation 
and explicitly point our attention on the physical assumptions disregarded elsewhere. We 
show that the usual assumptions can be relaxed leading to generalized formulas that hold 
even in the case of horizontal winds. This fundamental physics has some relevance to the 
current discussions on the climate debate. 
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1 Introduction 



In the following, we derive approximate temperature profiles of an atmosphere, also called 
adiabatic lapse rates or better barometric formulas. Our variant of a derivation goes beyond 
the common standard exercise of a thermodynamics lecture, where commonly the discussion 
of the underlying physical assumptions is missed. We depart from the Navier-Stokes equation 
and explicitly point our attention on the physical assumptions disregarded elsewhere. By the 
way, this derivation is a good example on how to apply the magnetohydrodynamic equations 
regarded as redundant by some of our critics. Furthermore, it explicitly shows that in physics 
an application of formulas is valid only in a finite space-time region. In addition, we show 
that the usual assumptions can be relaxed leading to generalized formulas that hold even in 
the case of horizontal winds. 

A brief historical review of the barometric formula is given in Ref. [1]. The reader is also 
referred to the textbook by Riegel and Bridger on "Fundamentals of Atmospheric Dynamics 
and Thermodynamics" [2]. 

2 On the derivation of the barometric formulas 
2.1 Input from hydrodynamics 

As described in our falsification paper [4, 3] the core of a climate model must be a set 
of equations describing the equations of fluid flow, namely the generalized Navier-Stokes 
equations. They describe the conservation of momentum and read 



where v is the velocity vector field, p the pressure field, $ the gravitational potential, R the 
friction tensor, and F cxt are the external force densities, which could describe the Coriolis 
and centrifugal accelerations. Neglecting the friction term and the electromagnetic fields we 
obtain the Euler equations. 

Assumption 1 

• We neglect the electromagnetic field terms. 

We get the more common version of the Navier-Stokes equations 



d_ 
di 



(gv) + V- (£v® v) = -Vp- gV<5> + QeE+j x B + V - R + F cxt 



(1) 



d_ 
dt 



(gv) + V • (^v® v) = -Vp-£V$ + V-R + F cxt 



(2) 
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The left hand side of this equation may be rewritten according to 
d dp <9v 

— (g v) + V • (q v <g> v) = -^v + g ■ — + V • (gv) v + gv ■ V v (3) 

With the continuity equation for the mass density dg/dt + V • (gv) = this term simplifies 
to 

£-^T + £v-Vv (4) 

Thus we obtain the well-know form of the Navier-Stokes equations, or, preferring the singular 
form, the Navier-Stokes equation 

Q- -oT + ffv Vv = -Vp-^V$ + V-R+F ext (5) 

at 

where the term — g V$ is gravity. If we neglect the viscosity term, we are left with the Euler 
equation 

<9v 

£-^7 + £v- Vv = -Vp-£V$ + F cxt (6) 
at 

Assumption 2 

• We assume that v(r, t) is independent of r. In sharp distinction to the standard deriva- 
tion of the barometric formulas, we relax the usual condition v(r, t) = in order to 
allow non- vanishing velocity fields v(r, t), which are independent of r. 

Consequently, the viscosity tensor R and the non-linear term f?v ■ Vv are zero, such that 

Q--q^ = - V P - ^ V$ + F cxt (7) 

Remark: If one writes 

Q v ■ Vv = -gv x (V x v) + q V Ql v | 2 ) ( 8 ) 

one could weaken this assumption to potential velocity fields. With these formulas one can 
derive the Bernoulli equation. 

In case of a rotating atmosphere of the Earth the last term F ex t of the right hand side of 
Eq. 7 describes the centrifugal acceleration and the Coriolis acceleration. The latter vanishes 
for a identically vanishing velocity field. 
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Assumption 3 

• We set F cxt to zero. 

We now have two fields, namely — Vp and — g V$, which will accelerate the volume elements, 
if they are different fields: 



Q 



dv 
~dt 



-Vp - g V$ 



Let us follow the common notation and write for the gravitational field 

g = -V$ 

Assumption 4 

• We assume that, as usual, acceleration due to gravity is vertical, i.e. we set 

g = -9 e z 



(9) 



(10) 



(11) 



where e z is the unit vector in z-direction and g is constant in space and time. This is 
the flat earth hypothesis. Furthermore, we neglect the variation of the gravitation field 
induced by the gravitation fields of the Sun, the Moon, and the planets. 

Assumption 5 

• We assume that the wind blows only horizontally, i.e. 



v • e z = 



(12) 



Eq. 9 now becomes 



^ g(dv x /dt) \ ( dp/dx ^ 



g (dvy/dt) 




dp/dy 
y dp/dz J 



( ° ^ 



V eg ) 



(13) 



V u / V dp l dz ) \ Qg ) 

That is, with the usual assumptions about geometry we would get the hydrostatic equation 



dp 

dz 



(14) 



without the usual assumption v(r,t) = 0. In standard thermodynamics for a macroscopic 
volume the pressure p is characterized by one number, not a field. Irreversible thermodynamics 
is a (classical) field theory and hydrodynamics is a special case. 
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2.2 Input from thermodynamics 

The equation of state for an ideal gas reads 

pv = RT (15) 

where v is the volume of one gram, and R = R/(l Mol), where R is the usual molar gas 
constant. Using the density g we also may write, respectively, 

I = RT, g=J- (16) 
Q RT 

Assumption 6 

• The air of the atmosphere obeys an equation of state of an ideal gas. 

With q = p/(RT) inserted into Eq. 14 we have 

dp = _pg_ = (lMol)pg 

dz RT ' RT 
If the molecular mass of the gas is greater, then the decrease of pressure with increase of 

height will be greater as well. For a temperature field T that is constant in space and time 

this equation can be integrated. 

2.3 The isothermal atmosphere 
Assumption 7a 

• We postulate an isothermal atmosphere. 
Separation of variable gives 



which may be integrated to 
yielding 

v = On exD I - 

RT 

from which, with g — p/ {RT), one obtains the density as a function of height 

6 RT * V RT s "V V RT' 

Thus, with help of these relations and assumptions, we obtain the barometric height formulas 

in case of an isothermal atmosphere. The lapse rates for the pressure and density, respectively, 

depend on the molecular mass of the gas, since R = R/(l Mol). 



p = p exp ( —£-(z - zq) ) (20) 



exp ("ir^ - ^) = ^° exp ("ir^ - ^) ^ 
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2.4 The adiabatic atmosphere 

In what follows, we need three relations for the heat differential form dQ, namely 

RT 

dQ = CJT)dT+ dv (22) 

v 

~ RT 

dQ = (C v (T) + R)dT dp (23) 

V 

dQ = ^(c v (^+R}pdv + lc v (^)vdp (24) 

Assumption 7b 

• For the ideal gas we use the reversible work form pdv. 

We now calculate die adiabatic state changes, i.e. we set dQ = 0, and, separating the variables, 
we obtain 

T\ til J vi V W2/ 

T 'm±* d T = f * = !„(£») (26) 
'n rt J pi p \Pi) 

Assumption 7c 

• The specific heats of an ideal gas are independent of the absolute temperature. 
Thus, we can continue our calculations with letters 



RJtiT r vrj \v 2 

C v + R f r,dT C. + 5 ln /^ =]n /M j 



R M T R \TJ \ Pl 

For constant heat capacities we integrate the third equation (24) for dQ = 0. Since 

= l?(C v + R)pdv + i*C v vdp (29) 
R R 

is equivalent to 

= (C v + R)p dv + C v v dp (30) 

we get 

~ f V2 dv f'P 2 dv 

(C v + R) — = -C v ^ (31) 

J vi V J P1 p 
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C " + R In (*) =j!L\ (32) 



and, therefore, 

C 

C v '"\viJ '"\P2, 

With our assumptions, we have 

Cp — C v = R (33) 

Setting 

• C p /C v = K, 

• T 2 — T,T 1 — T , 

• P2=P,Pl= PO, 

• v 2 = v,v 1 = v , 

we obtain the well-known adiabatic equations of state: 

K .h« = ln(V| (34) 



v Q j \ p 

= * (35) 

p / V 

Po \vo 



and 





= In | 


fv 






V v 




= In | 








V v 



K-l 

Ty/(-D ^ 
T J v 



and 



K-l VT / \po 



P 



rp > k/(k—1) 



(36) 



C V-M^r) = In(-) (37) 



' M^r) = M~) (38) 



(39) 



* =m(l] (41) 



(42) 



These adiabatic equations of state are well-known from standard textbooks. However, almost 
never the assumptions are discussed, under which they hold true. 
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p_ _ 




Po 




Q_ 




Qo 


(I 


p_ 


(I 


Po 





Now we replace the one-gram volume v by the density p. We get 

(43) 

(K-l) 

(44) 

(k-1) 

(45) 

Rewriting the first of these three equations as 

'-(I)' (46) 

and inserting it into Eq. 14, namely dp/dz = —gg, whereby we consider p and g as functions 
of z, we obtain, applying the chain rule of differential calculus, 

dp Po K -idg , * 

^ = V Tz = (47) 

= -9 (48) 

Qq dz 

±{(T X ) = - {K - l)90 ° (49) 
dz k po 

Integrating this equation we obtain the density g as a function of height 

g = g [z-zo) (50) 

upo 

Q{z) = , fl-^lM(,_, o) y /M (51) 

From this and Eq. 46 we get the pressure p as a function of height z 

(i i \ \ «/(«— i) 

1 - ^lM (2 _ Zo) ( 52) 

Inserting this into p/p = (T/T ) K ^ K ~^ resp. T = T (p / po)^ 1 ^ K we obtain the temperature 
T as a function of the height z 

T{z)=T 1 (z-^o) = T o ~ — (2-20) (53) 

V K Po J kR 

We conclude: 



• The temperature is decreasing linearly with increasing height. 
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As in the case of an isothermal atmosphere the lapse rate will be higher if the molecular mass is 
greater. In the adiabatic case this slope depends on the adiabatic coefficient k. Note, that one 
can always weaken the assumption to get similar results (e.g. moist adiabatic approximation 
and so on). 

With our assumptions we may rewrite the constants appearing in the formula before the 
height term, i.e. 



(k - l)g g (C p /C v -l)g g 
Kpo (C p /C v )p 
(C p - C v )g g 

CpPo 
(C p -C v )g 

CpRT 
9 



and 



C p T 

[k - 1) g _ (Cp -C v )g _ g 



kR CpR Cp 

Finally, with these coefficients the adiabatic height formulas read 



(54) 

(55) 

(56) 
(57) 

(58) 



/ \ Cv j ( Cp Cv ) 

g(z) = gJl--^-( z - ZQ )\ (59) 



C p Tq 
9 

C p Tq 



/ \ Cp/(Cp—Cv) 

p(z) = p (l--J-( z - Zo )) (60) 



T(z) = T -Mz-zo) (61) 
Op 

3 Results 

By combining hydrodynamics, thermodynamics, and imposing the above listed assumptions 
for planetary atmospheres one can compute the temperature profiles of idealized atmospheres. 
In case of the adiabatic atmosphere the decrease of the temperature with height is described 
by a linear function with slope —g/C p , where C p depends weakly on the molecular mass. As 
elucidated in our paper [4, 3] mixtures of gases are covered in the context of Gibbs ther- 
modynamics. Since the measurable thermodynamic quantities of a voluminous medium, in 
particular the specific heat and the thermodynamic transport coefficients, naturally include 
the contribution from radiative interactions, we cannot expect that a change of concentration 
of a trace gas has any measurable effect. At this point, it is important to remember that the 
barometric formulas do not hold globally but have only a limited range of validity. 
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Appendix: Relevance to the current climate debate 

In our falsification paper [4, 3] we have shown that the atmospheric C0 2 greenhouse effects 
[5] as taken-for-granted concepts in global climatology do not fit into the scientific framework 
of theoretical and applied physics. By showing that 

(a) there are no common physical laws between the warming phenomenon in glass houses 
and the fictitious atmospheric greenhouse effects 

(b) there are no calculations to determine an average surface temperature of a planet 

(c) the frequently mentioned difference of 33 degrees Celsius is a meaningless number cal- 
culated wrongly 

(d) the formulas of cavity radiation are used inappropriately 
(f) the assumption of a radiative balance is unphysical 

(e) thermal conductivity and friction must not be set to zero 

the atmospheric C0 2 greenhouse effects have been refuted within the frame of physics [4, 3]. 

In other words, the greenhouse models are all based on simplistic pictures of radiative 
transfer and their obscure relation to thermodynamics, disregarding the other forms of heat 
transfer such as thermal conductivity, convection, latent heat exchange et cetera. Some of 
these simplistic descriptions define a "Perpetuum Mobile Of The 2nd Kind" and are therefore 
inadmissable as a physical concept. 

In the speculative discussion around the existence of an atmospheric natural greenhouse 
effect [6] or the existence of an atmospheric C0 2 greenhouse effect it is sometimes stated that 
the greenhouse effect could modify the temperature profile of the Earth's atmosphere. This 
conjecture is related to another popular but incorrect idea communicated by some proponents 
of the global warming hypothesis, namely the hypothesis that the temperatures of the Venus 
are due to a greenhouse effect. For instance, in their book "Der Klimawandel. Diagnose, 
Prognose, Therapie" (Climate Change. Diagnosis, Prognosis, Therapy) "two leading inter- 
national experts", Hans- Joachim Schellnhuber and Stefan Rahmstorf, present a "compact 
and understandable review" of "climate change" to the general public [8]. On page 32 they 
explicitly refer to the "power" of the "greenhouse effect" on the Venus. 

The claim of Rahmstorf and Schellhuber is that the high venusian surface temperatures 
somewhere between 400 and 500 Celsius degrees are due to an atmospheric C0 2 greenhouse 
effect [8]. Of course, they are not. On the one hand, since the venusian atmosphere is 
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opaque to visible light, the central assumption of the greenhouse hypotheses is not obeyed. 
On the other hand, if one compares the temperature and pressure profiles of Venus and 
Earth, one immediately will see that they are both very similar. An important difference is 
the atmospheric pressure on the ground, which is approximately two orders higher than on 
the Earth. At 50 km altitude the venusian atmospheric pressure corresponds to the normal 
pressure on the Earth with temperatures at approximately 37 Celsius degrees. However, 
things are extremely complex (volcanic activities, clouds of sulfuric acid), such that we do 
not go in details here [9]. 
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